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1 Introduction 


In general relativity and cosmology, our knowledge about spatially homogeneous cosmo¬ 
logical models has increased substantially over the years, and we are able to say that, for 
a large number of models, the qualitative behaviour of solutions is now well understood; 
see []P for an overview. The majority of results, however, concerns solutions of the Einstein 
equations coupled to a perfect fluid (usually with a linear equation of state). It is thus 
important to note that these results are in general not robust, i.e., not structurally stable, 
under a change of the matter model; significant changes of the qualitative behaviour of 
solutions occur, for instance, for collisionless matter. 

Several fundamental results on spatially homogeneous diagonal models of Bianchi type I 
with collisionless matter have been obtained in |2j. Diagonal locally rotationally symmetric 
(LRS) models have been investigated successfully by using dynamical systems methods, 
see ini for the case of massless particles, and HE] for the massive case. In particular, 
solutions have been found whose qualitative behaviour is different from that of any perfect 
fluid model of the same Bianchi type. 

The purpose of this article is to re-investigate the diagonal (non-LRS) Bianchi type I 
models with collisionless matter. Our analysis is based on dynamical systems techniques, 
which enables us to obtain a more detailed picture of the global dynamics than the one 
previously given. In particular, we show that the dynamical behaviour toward the past 
singularity of the collisionless matter model differs considerably from the Bianchi type I 
perfect fluid model. 

The outline of the paper is as follows. In Section[2]we reformulate Einstein’s field equations 
for the diagonal Bianchi type I case with collisionless matter as a reduced dimensionless 
dynamical system on a compact state space. In Section O) we give the fixed points of the 
system and list and discuss a hierarchy of invariant subsets of the state space, which is 
associated with a hierarchy of monotone functions. In SectionHwe hrst present the results 
of the local dynamical systems analysis; subsequently, we focus on the analysis of the 
global dynamics: we establish two theorems that formulate the future and past asymptotic 
behaviour of solutions, respectively. As regards the future asymptotics we show that 
all models isotropize asymptotically toward the future. The past asymptotic behaviour 
is more complicated since there exists several types of past asymptotic behaviour; in 
particular we establish that the past attractive set resides on a set that contains a so- 
called heteroclinic network. The proof of the theorems is based on methods from global 
dynamical systems analysis; in particular, we exploit the hierarchy of monotone functions 
in conjunction with the monotonicity principle. Finally, in Section E) we conclude with 
some remarks about our results and their implications. Appendix HI provides a brief 
introduction to relevant background material from the theory of dynamical systems; in 
particular we cite the monotonicity principle. The proofs of some of the statements in 
the main text are given in Appendix H and o In Appendix El we discuss the physical 
interpretation of one of the most important boundaries of our state space formulation. 
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2 The reflection-symmetric Bianchi type I Einstein-Vlasov 
system 

In a spacetime with Bianchi type I symmetry the spacetime metric can be written as 

ds'^ = —dt^ + gij{t)dx^dx^ , i,j = 1,2,3 , (1) 

where gij is the induced Riemannian metric on the spatially homogeneous surfaces t = 
const. Since the metric is constant on t = const, it follows that the Ricci tensor of g^j 
vanishes. Einstein’s equations, in units G = 1 = c, decompose into the evolution equations. 


dtgij = -2kij , dtk^j = irk feE - SnT) + - p) - , (2a) 

and the Hamiltonian and momentum constraint 

(tr/c)^ — k^jk\ — IOttp — 2A = 0 , jk = 0 ■ (2b) 

Here, kij denotes the second fundamental form of the surfaces t = const. The matter 
variables are defined as components of the energy-momentum tensor T^i, (p, = 0,1,2,3), 


according to p = Tqo, jk = Tofe! denotes the spatial components. The cosmological 
constant A is set to zero in the following; the treatment of the case A > 0 is straightforward 
once the case A = 0 has been solved, cf. the remarks in the conclusions. 

In this paper we consider collisionless matter (Vlasov matter), i.e., an ensemble of freely 
moving particles described by a non-negative distribution function / defined on the mass 
shell bundle PM C TM of the spacetime; for simplicity we consider particles with equal 
mass m. The spacetime coordinates (t, x*) and the spatial components u* of the four- 
momentum (measured w.r.t. d/dx^) provide local coordinates on PM so that / = 
f{t,x^,v^). Compatibility with Bianchi type I symmetry forces the distribution function 
/ to be homogeneous, i.e., / = f{t, v^). The evolution equation for / is the Vlasov equation 
(the Liouville equation) 

dtf + = dtf + 2k\v^d,,! = 0 . (2c) 

The energy-momentum tensor associated with the distribution function / is given by 

= I fvf^v^'volpM , 

where volpM = {det g)^/‘^VQ^dv^dv^dv^ is the induced volume form on the mass shell; 
vq is understood as a function of the spatial components, i.e., Uq = -|- gijv^vA The 

components p, jk, and Tij, which enter in (1^^ and (I2bl) can thus be written as 


p = ^ 

^ / {m^ + {deig)~^/‘^dvidv2dv^ , 

(2d) 

jk 

^ fvk{deig)~^/‘^dvidv2dv2, , 

(2e) 

“ j 

1" fviVj (m^ + g’^^VkViJ ^ {d(i.ig)~^l‘^dv\_dv2dvi . 

(2f) 
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The Einstein-Vlasov system m is usually considered for particles of mass m > 0, however, 
the system also describes massless particles if we set m = 0. (For a detailed introduction 
to the Einstein-Vlasov system we refer to |Hj and (Hj.) 

The general spatially homogeneous solution of the Vlasov equation (123) in Bianchi type I 
is 

= (3) 

where the vi are the covariant components of the momenta and /o is an arbitrary non¬ 
negative function, see HD]. (By inserting (|3]) into (1^ and using that Vi = gij{t)v^ it is 
easy to check that fo{vi) is a solution.) The momentum constraint in pb]) then reads 


j fo(vi)vkdvidv2dv3 = 0 . 


( 4 ) 


Henceforth, for simplicity, /o is assumed to be compactly supported. 

There exists a subclass of Bianchi type I Einstein-Vlasov models that is naturally associ¬ 
ated with the constraint (HJ: the class of “reflection-symmetric” (or “diagonal”) models. 
The following symmetry conditions are imposed on the initial data: 

fo{vi,V2,V3) = fo{-Vi,-V2,V3) = fo{-Vl,V2,-V3) = fo{vi,-V2,-V3) , (5a) 

gij{to) ,kij{to) diagonal. (5b) 

These conditions ensure that Tjj(fo) is diagonal, hence gij, kij, and Tij are diagonal for 
all times by ©• In the present paper, we will be concerned with this class of reflection- 
symmetric models. 

The Einstein-Vlasov system m thus reduces to a system for six unknowns, the diagonal 
components of the metric gii{t) and the second fundamental form fe*j(t) (no summation). 
The equations are (123 and the Hamiltonian constraint in (|2E|). The initial data consists 
of gii{to), /c\(to); in addition we prescribe a distribution function fo{vi) that provides the 
source terms in the equations via (|2H|) and m- 

In the following we reformulate the Einstein-Vlasov system as a dimensionless system on 
a compact state space. To that end we introduce new variables and modified matter 
quantities. Let 

(< 5 ) 

i 

and define the dimensionless variables 


g^i 

hi 

Si :=-^ - 1 , 

H 

9 

m 

(7a) 

Si ■ — ) 

X 

' 9 1 ’ 

+ X 

where s* = 1 , 

i 

M 

II 

o 


(7b) 

The transformation from the variables 
to the above constraints, is one-to-one. 

{gii,k\) to {si,T,i,x,H), where (si,Si) are 
(Note that x can be obtained from 2; when 

subject 
m > 0.) 
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By distinguishing one direction (1, 2, or 3), one can decompose the Si and simultaneously 
introduce a trace-free adaption of the shear to new Sj- variables as is done in, e.g., PP; 
however, since Bianchi type I does not have a preferred direction we will refrain from doing 
so here. 


Next we replace the matter quantities p, T\ (no summation) by the dimensionless quan¬ 
tities 


n : = 


Sirp 



w : = 



^EiT\ 
3 P 


( 8 ) 


Expressed in the new variables, Wi can be written as 


Wi = 


(1 - z)si j fo vf [z + {1- z) Sk vl] dvidv 2 dv 3 
//o [^ + (1 - z) Sk vl] dvidv 2 dv 3 


(9) 


Finally, let us introduce a new dimensionless time variable r defined by 


dr = H ^dt ; 


( 10 ) 


henceforth a prime denotes differentiation w.r.t. r. 

We now rewrite the Einstein-Vlasov equations as a set of dimensional equations that 
decouple for dimensional reasons and a reduced system of dimensionless coupled equations 
on a compact state space. The decoupled dimensional equations are 


n 


1 —(1 -" w ) 
2 ^ ' 


H' = -3H 
x' = -2x ( 1 -f 


(lla) 

(llb) 


The reduced dimensionless system consists of the Hamiltonian constraint, cf. (f2E|l . 


1 - H = 0 , where := 1 ^ , 

k 

which we use to solve for H, and a coupled system of evolution equations 


-311 

^(1 - w)Y.i - 

[Wi - w) 

-2si 




k 


z = 2z {I - z) i 1 -F Sfc Efc 

V k 



( 12 ) 


(13a) 

(13b) 

(13c) 


In the massive case m > 0 the decoupled equation for x is redundant since the equation for 
z is equivalent. In the massless case m = 0 we have z = 0; hence, although the equation 
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for X does not contribute to the dynamics, x is needed in order to reconstruct the spatial 
metric from the new variables. 

The dimensionless dynamical system m together with the constraint m describes the 
full dynamics of the Einstein-Vlasov system of Bianchi type I. In the massive case the 
state space associated with this system is the space of the variables {(Sj, Sj, z)}, i.e., 

T := {(Ei, Si, z) I < l) A (si > 0) A (0 < 2 ; < 1)} , (14) 

where the Si and Ej are subject to the constraints = 1, Yl,k = 0. (The inequal¬ 

ities for Si and Ej follow from the definition (I7al) and the constraint (1121) . respectively.) 
The state space X is thus five-dimensional. 

It will turn out eventually that all solutions asymptotically approach the boundaries of 
X-. z = 0, z = 1, Si = 0, Q = 0 (<tA T? = 1). This suggests to include these sets in the 
analysis, whereby we obtain a compact state space X. 

The equations on the invariant subset z = 0 of X are identical to the coupled dimensionless 
system in the case of massless particles m = 0. We will therefore refer to the subset z = 0 
as the massless subset; it represents the four-dimensional state space for the massless case. 

We conclude this section by looking at some variables in more detail. The inequality E^ < 1 
together with the constraint = 0 results in |Ej| < 2 for all i. Note that equality 

is achieved when (Ei,E 2 ,E 3 ) = (± 2 ,:pl, 2 pl) and permutations thereof, cf. Figure^ The 
matter quantities satisfy 

0<tc<|, 0 < Wi < 3w < 1 . (15) 

The equalities hold at the boundaries of the state space: tc = 0 iff 2 ; = 1, and w = ^ 
iff 2 ; = 0 ; tCj = 0 iff Si = 0 , and Wi = Stc iff s* = 1 (provided that 2 ; < 1 ; for 2 ; = 1 , 
0 = rcj = Stc = 0). 

There exists a number of useful auxiliary equations that complement the system (USD: 


n' = n 


p =-p 


3(1 - w)T.‘^ -'^Wk^k 


[3{l + w) x'^Wk^k] < -2/>. 


(16) 

(17) 


The inequality in CZl) follows by using Ej > — 2 Vi and m- This shows that p increases 
monotonically toward the past which yields a matter singularity, i.e., p —> 00 for r —> 
— 00 . It is often beneficial to consider the equations of the original variables as auxiliary 
equations, e.g., ( 5 **)' = — 25 :®* (1 -|- Ej). 


3 Fixed points, invariant subsets, and monotone functions 

3.1 Fixed points 

The dynamical system m possesses a number of fixed points, all residing on the bound¬ 
aries of the state space; see Tabled 
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Fixed point set 

defined by 

FS^ 

z = l,^j = 0 Mj 

KCi 

z = 1, T? = Si = 1 , Sj = 0 Vj 7 ^ i 

TSi 

0 < z < 1, T,i = 2, T,j = —1 Vj ^ i, Si = 0 

FU 

z = 0, 'Ej = 0 Vj, Wj = 1/3 Vj 


z = 0, Si = 0, Tii = — 1 , Sj = 1/2 = Wj Vj 7 ^ i 

QL? 

z = 0, T,i = — 2 , Sj = 1 Vj 7 ^ z. Si = 0 

KC» 

z = 0 , S^ = 1 , Sj = 1 , Sj = 0 Vj 7 ^ z 


Table 1: The fixed point sets. The range of the index i is always i = 1... 3. The superscript 
denotes the value of z; the first kernel letter describes the type of fixed point set; if there 
is no second kernel letter the fixed point set is just a point; if there is a second kernel 
letter this letter denotes the dimensionality and character of the set — S refers to surface, 
L stands for line, and C for circle. 


FS^ is a surface of fixed points that correspond to the flat isotropic dust solution. The 
circles consist of fixed points with constant Sj, satisfying = 0 , 

see Figure ^ these fixed points correspond to Kasner solutions. The fixed points on 
TSj are associated with the Taub representation of the flat Minkowski spacetime. The 
intersection of TSj with [z = 0) yields a line of fixed points which we denote TL?. The 
fixed points on QL° correspond to the non-flat LRS Kasner solutions. F*^ is a fixed point 
that corresponds to the flat isotropic radiation solution. In Appendix H we prove that 
F® is well-defined through the equations wi = W 2 = = !/?> (which are to be solved 

for (si, S 2 , S 3 )). The location of F° depends on the chosen distribution function, since the 
equations tci = rc 2 = tea = 1/3 involve /q. Analogously, the equations rcj = 1/2 (V j 7 ^ i) 
yield a unique solution, the fixed point D/; the location of the point D? also depends on 
/q. The fixed points D? are associated with a scale-invariant LRS solution (related to a 
distributional /o; see Appendix m for details). 

The LRS points on KC/ play a particularly important role in the following, which motivates 
that they are given individual names. We denote the three Taub points on KC^ defined 
by Sj = 2 (and thus S; = — 1 V/ 7^ j) by T/^, while we denote the three non-flat LRS 
point on KC/ given by Sj = —2 (and thus S; = 1 V/ 7^ j) by Q/^. The Kasner circles 
KC^ and KC° are connected by the lines TL/ and QL/; the end points of the line TL/ 
are the Taub points T^j and T°^; analogously, the end points of QL/ are the points 
and Q°j. (Here, {i,j, k) is an arbitrary permutation of (1, 2, 3).) The remaining points T|) 
and Q[) (/ = 1... 3) do not lie on any of the fixed point sets TL/ or QL/. This fixed point 
structure is depicted in Figure |2J 

3.2 Invariant subsets and monotone functions 

The dynamical system m possesses a hierarchy of invariant subsets and monotone func¬ 
tions. Since this feature of the dynamical system will turn out to be of crucial importance 
in the analysis of the global dynamics, we give a detailed discussion. 
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JY: On the full (interior) state space X define 

M(i) = • ( 18 a) 

A straightforward computation shows 

M(\) = 2 M(i) , (18b) 

i.e., M(i) is strictly monotonically increasing along orbits in X. Note that is inti¬ 
mately related to the spatial volume density since det( 5 fjj)^/^. 

Z^: This subset is characterized by z = 1. Since Wi = w = 0, the equations for Sj decouple, 
and the essential dynamics is described by the equations S' = — (3/2)(l — S^)Sj. (Note 
that these equations are identical to the Bianchi type I equations for dust — it is therefore 
natural to refer to Z^ as the dust subset.) Explicit solutions for these equations can be 
obtained by noting that Si oc S 2 oc S 3 for all solutions, or by using that fl' = 3S^ Q. 

Z^: This subset is the massless boundary set z = 0. Since w = 1/3, the dynamical 
system m reduces to 


S' = -11 [ 1 -h Si - 3t(;i ] , s' = -2si [ Si - Sk S 
Consider the function 

,lV2 


(19) 


M( 2 ) = (1 - S^) ^ (S1S2S3) I fo 1^^^ Skvl 


dvidv2dvs 


( 20 a) 







Figure 2; A schematic depiction of the fixed points on z = 0. The underlying structure is 
the three sides of the Sj-triangle si + S 2 + = 1: each point represents a disc < 1; the 

vertices contain the Kasner circles KC^. Bold lines denote the lines of fixed points TL®, 
QL°, and KC^ 


The derivative is 

M(' 2 ) = - 2 S 2 M( 2 ) , ( 20 b) 

which yields monotonicity when ^ 0. If = 0, then 

M' 2 )= 0 , M'')=0, M'" =- 6 M( 2 ) (20c) 

i 

Hence, M( 2 ) is strictly monotonically decreasing everywhere on z = 0, except at the fixed 
point (for which = 0 and wi = W 2 = W 3 = 1/3), where M( 2 ) attains a (positive) 
minimum. The latter follows from the fact that (1 —S^)“^ is minimal at the point Sj = OVi 
and that dM(^ 2 )/dsi = { 2 si)~^[wi — 1 / 3 ]M( 2 ). 

Si {i = 1,2,3): These invariant boundary subsets are defined by s* = 0 (which yields 
Wi = 0). There exists a monotone function on 5i, 

M(3) = (s2S3)-'/2 , M('3 ) = (2 - Si) M(3) ; (21) 

analogous functions can be obtained on S 2 and ^3 through permutations. 

JC: This boundary subset is the vacuum subset defined by H = 0 (or equivalently S^ = 1). 
The Sj are constant on this subset, which completely determines the dynamics of the Si 
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variables (via JIsBl) or via the auxiliary equation for The Bianchi type I vacuum 

solution is the familiar Kasner solution and we thus refer to /C as the Kasner subset. 

Intersections of the above boundary subsets yield boundary subsets of lower dimensions; 
those that are relevant for the global dynamics are discussed in the following. 

5? and S}: The intersection between the subset Si and and Z^ yields three-dimensional 
invariant subsets (sj = 0) H (z = 0) and (sj = 0) H (z = 1), respectively. On there exists 
a monotonically decreasing function: 

= (1 + , -^( 4 ) ~ . ( 22 ) 


Sij-. These subsets are defined by setting Sj = 0 and Sj = 0 (j 7 ^ i), i.e., Sij =5^0 Sj. On 
Sij, we have = 1 (A: 7 ^ i,j) and Wk = 3w, because Wi = Wj = 0. 

T)^: The subsets 5? admit two-dimensional invariant subsets characterized by {z = 
0) n (sj = 0) n (Bj = — 1). On consider the function 

M( 5 ) = {2 + 1 : 21 : 3 )“^ ( 5253 )“^^^ J fo [S 2 V 2 + dvidv 2 dv 3 ; (23a) 

analogous functions can be defined on and Eqs. (jlDj) imply 





(5) 


(1-2B2)V(1-2B3)' , 


(23b) 


i.e., is strictly monotonically decreasing unless B 2 = 1/2 = B 3 . In the special case 
B 2 = 1/2 = B 3 we obtain 


V,'5)=0. “",1=0. “l"l = -¥“( 


'(5) 


'(5) 


(5) 


{W2 - + {W3 - 


(23c) 


Hence, is strictly monotonically decreasing everywhere on except for at the fixed 
point Di, for which T ,2 = '^3 = W 2 = W 3 = cf. Tabled The function M( 5 ) possesses a 
positive minimum at Di. This is because ( 2 -|-I] 2 B 3 )“^ is minimal at the point B 2 = B 3 = 
1/2 and dMf^^'^/dsi = (2sj)“^[r(;i — 1 / 2 ]M( 5 ) for i = 2,3. 

KP\ The intersection of the Kasner subset K, = (B^ = 1) with the z = 0 subset yields a 
3-dimensional subset, /C®. This subset will play a prominent role in the analysis of the 
past asymptotic behaviour of solution. 

The remaining subsets we consider are not located at the boundaries of the state space, 
but in the interior; these subsets are invariant under the flow of the dynamical system, if 
the distribution function /o satisfies certain symmetry conditions. 

LRSi: We define the subset LRSi of X through the equations B 2 = B 3 , ^2 = u) 3 ; LRS 2,3 
are defined analogously. In order for these sets to be invariant under the flow of the 
dynamical system, the distribution function /o must satisfy conditions that ensure com¬ 
patibility with the LRS symmetry, see Appendix^for details. For an orbit lying on LRSi, 
Equation ()13b|l entails that S2{t) oc S3(r) (where the proportionality constant exhibits a 
dependence on /o, which enters through the equation W 2 = W 3 ), and hence §22 oc 533 ; by 
rescaling the coordinates one can achieve 522 = 533 , he., a line element in an explicit LRS 


10 






form. Hence, the LRSj subsets, if present as invariant subsets, comprise the solutions with 
LRS geometry. 

FRW: If /o is compatible with an isotropic geometry, see Appendix m for details, the 
one-dimensional subset characterized by the equations Sj = 0 Vi and wi = W2 = W3 = w 
is an invariant subset (in fact: orbit), the FRW subset. The equations Sj = 0 yield 
Si = const, whereby we obtain a Friedmann-Robertson-Walker (FRW) geometry, since the 
spatial coordinates can be rescaled so that gtj oc dij. Note that the location in Si of the 
FRW subset depends on /o, since the equations wi = W 2 = 103 , which are to be solved for 
(si, 52 , 53 ), involve /o, cf. AppendixEl Remarkably, in the massless case the existence of 
a FRW solution (which corresponds to the fixed point F°) does not require any symmetry 
conditions on /q. 


4 Local and global dynamics 

4.1 Local dynamics 

Let us consider smooth reflection-symmetric Bianchi type I Vlasov solutions that approach 
fixed point sets when r ^ — 00 . 

Theorem 4.1. In the massive (massless) case there exists 

(a) a single orbit that approaches (corresponds to) F°, 

(b) three equivalent one-parameter sets of orbits (three single orbits) that approach D?, 
f = 1...3, 

(c) one three-parameter (two-parameter) set of orbits that approaches QL 5 ; QL® and 
QL 3 yield equivalent sets, 

(d) one four-parameter (three-parameter) set of orbits that approaches the part o/KC 5 
defined by 1 < T,i < 2 ; KC 2 and KCg yield equivalent sets. 

Proof. The statements of the theorem follow from the local stability analysis of the fixed 
point sets F^, D?, QL?, KC^, when combined with the Hartman-Grobman and the reduc¬ 
tion theorem, since the fixed points F*^, D? are hyperbolic and QL®, KC^ are transversally 
hyperbolic. This requires the dynamical system to be and this leads to some restric¬ 
tions on /q. However, it is possible to obtain an alternative proof that does not require 
such restrictions. Such a proof can be obtained from the hierarchical structure of invariant 
sets; we will refrain from making the details explicit here, since our analysis of the global 
dynamics below contains all essential ingredients implicitly. 

Interpretation of Theorem EH (massive case): A three-parameter set of solutions con¬ 
verges to every non-LRS Kasner solution as t —> 0. (In the state space description three 
equivalent sets of orbits approach three equivalent transversally stable Kasner arcs that 
cover all non-LRS Kasner solutions; the equivalence reflects the freedom of permuting 
the coordinates.) Furthermore, a three-parameter set of solutions approaches the non-flat 
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LRS Kasner solution. Hence, in total, a four-parameter set of solutions asymptotically ap¬ 
proaches non-flat Kasner states. There exist special solutions with non-Kasner behaviour 
toward the singularity: one solution isotropizes toward the singularity and a one-parameter 
set of solutions approaches a non-Kasner LRS solution of the type (three equivalent 
one-parameter sets of orbits approach three equivalent non-Kasner LRS fixed points as¬ 
sociated with this solution). For the latter solutions = 3/4; these solutions cannot be 
interpreted as a perfect fluid solutions since they possess anisotropic pressures. 

In the following we show that the list of Theorem 14.1 1 is almost complete: there exist no 
other attractors toward the singularity with one exception, a heteroclinic network that 
connects the flat LRS-Kasner points. 

4.2 Global dynamics 

Theorem 4.2. All orbits in the interior of the state space X of massive particles [state 
space of massless particles] converge to FS^ [F^] when r ^ -|-oo; i.e., all smooth 
reflection-symmetric Bianchi type I Vlasov solutions isotropize toward the future. 

A proof of this theorem has been given in j^. In Appendix [Owe present an alternative 
proof based on dynamical systems techniques. 

Theorem 4.3. The a-limit set of an orbit in the interior of the state space is one of the 
fixed points F*^, D/, QL/, KC/, see Theorem \4.1\ or it is the heteroclinic network Bfi. The 
a-limit set of a generic orbit resides on the union of the fixed point sets KC/ and possibly 
the heteroclinic network Bfi. 

Remark. A heteroclinic network is defined as a compact connected flow-invariant set that is 
indecomposable (all points are connected by pseudo-orbits), has a finite nodal set (the set 
of recurrent points is a finite union of disjoint compact connected flow-invariant subsets), 
and has finite depth (taking the a/w-limit iteratively yields the set of recurrent points 
after a finite number of iterative steps); for details see and references therein. A 
simple case is a heteroclinic network of depth 1 whose nodes are equilibrium points: it can 
be regarded as a collection of entangled heteroclinic cycles. The heteroclinic network TC^ 
is of the latter type; it will be introduced in the proof of the theorem. 

The remainder of this section is concerned with the proof of Theorem 14.31 The first step 
in the proof is to gain a detailed understanding of the dynamics on the relevant invariant 
subspaces of the dynamical system. 


Dynamics on 5? 

Lemma 4.4. Consider an orbit in the interior of . Its a-limit set is either a fixed point 
on KCj or KC° (iff f k), QL? or TL/, or it is the heteroclinic cycle ifi]. The uj-limit 
set consists of the fixed point D/. 

Remark. The heteroclinic cycle will be defined in (EH). 
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Figure 3: Flow on the boundaries and on the invariant subset Si = —1 on 5i. The fixed 
point on Si = —1 is the point Dj*; the heteroclinic cycle TL^ consists of the fixed points 
T 22 , T 32 , T 33 , T 23 , and the connecting orbits. 


Proof. Without loss of generality we consider 5^, which can be described by the variables 

0 < S 2 < 1 (s 3 = 1 - S 2 ) and Si,S 2 ,S 3 . Sj = 0 , S^ < 1 ^ ; (24) 

hence 5i is represented by the interior of a cylinder, cf. Figure IHI The boundary of 5i 
consists of the lateral boundary 5° n the base 5 ^ 2 ) and the top surface 5 ^ 3 . 

Since 5i n is part of /C, it follows that Sj = const for all orbits on 5i n We 
observe that S 2 is monotonically increasing (decreasing) when S 2 < S 3 (S 2 > S 3 ), since 
S 2 = — 2 ^ 2(1 — S 2 )(S 2 — S 3 ); the two domains are separated by the lines of fixed points 
TL3 and QL3, see Figure El 

The key equations to understand the flow on 5 i 2 are 

n' = 0(2S2 - S 3 ) and S(j = 12(2 - S 3 ). (25) 

From the first equation it follows that all points on KC 3 are transversally hyperbolic 
repelling fixed points except for T 33 ; from the second equation we infer that T 33 is the 
attractor of the entire interior of 5 i 2 . Similarly, T 22 is the attractor on see Figure El 

The plane D)*, defined by Si = —1, is an invariant subset in 5i. In the interior of the 
plane we find the fixed point D^; the boundary consists of a heteroclinic cycle Tii, 

n-ZO . rpO rpO rpO rpO rpO 

'T-l ■ -*-22 ^ a32 ^ -‘-33 ^ -*-23 ^ -‘■22 • 

(Note that analogous cycles and exist on 5° and S^, respectively.) The function 
M( 5 ) is monotonically decreasing on T>^, cf. (1231) . Application of the monotonicity principle, 
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see Appendix^ yields that D? is the w-limit and that l-i\ is the a-limit for all orbits on 
Vl cf. Figure 13 

Consider now an orbit in 5]^ with Si 7 ^ —1. The function M( 4 ) = (1 + Si)^is monotonically 
decreasing on 5^, cf. ( 1221 ) . The monotonicity principle implies that the w-limit lies on 
Si = — 1 or S^ = 1 (but Si 7 ^ i 2 ). Since the logarithmic derivative of is positive 
everywhere on 5° n KP (except at T 22 and T 33 ), i.e., = 2 — '^^Wk^k > 0, it 

follows that the “wall” 5i H/C® of the cylinder is repelling everywhere away from Si = — 1 . 
Consequently, the w-limit of the orbit cannot lie on S^ = 1 but is contained in Si = — 1 . 
The fixed point is a hyperbolic sink, as we conclude from the dynamics on Si = — 1 
and from (1 + Si)“^(l + Si)'|qo = —3/4. Therefore, the a priori possible w-limit sets on 
Si = — 1 are Di and the heteroclinic cycle 1 -L\. 

To prove that Di is the actual tu-limit we again consider the function M( 5 ). However, we 
no longer restrict its domain of definition to but view it as a function on 5/; we obtain 

12M'5) = -M( 5 ) [(Si + 2S2)2 + (Si + + 6 (Si + 1)^ - 6 (Si + 1)] . (27) 

The bracket is positive when Si < —1; hence is decreasing when Si < —1. This 

prevents orbits with Si < — 1 from approaching TL^, since the cycle is characterized by 
M( 5 ) = 00 . Now suppose that there exist an orbit in Si > —1, whose w-limit is Ti'l. At late 
times the trajectory shadows the cycle; hence, for late times, the bracket in (EZI) is almost 
always positive along the trajectory - only when the trajectory passes through a small 
neighbourhood of (Si,S 2 ,S 3 ) = (— 1 , 1 / 2 , 1 / 2 ) the bracket is marginally negative. Since 
the trajectory spends large amounts of time near the fixed points and the passages from 
one hxed point to another become shorter and shorter in proportion, it follows that at late 
times is decreasing along the orbit (with ever shorter periods of limited increase). 
This is a contradiction to the assumption that the orbit is attracted by the heteroclinic 
cycle. We therefore draw the conclusion that Di is the global sink on 5i. 

Consider again an orbit in 5^ with Si 7 I —1. Invoking the monotonicity principle with the 
function M( 4 ) we find that the a-limit of the orbit must be located on S^ = 1 , Si 7 I — 1 . 
From the analysis of the flow on the boundaries of the cylinder we obtain that all fixed 
points on S^ = 1 except for T 22 and T 33 are transversally hyperbolic. The fixed points 
on KC2 with S2 < S3 and the points on KC3 with S2 > S3 are saddles; the fixed points 
on KC 2 with S2 > S3 and those on KC 3 with S2 < S3 are transversally hyperbolic 
sources (except for T 22 , T 33 ): every point attracts a one-parameter set of orbits from 5i 
as r —> — 00 . In contrast, each fixed point on TL 4 and QL 4 is a source for exactly one 
orbit. The structure of the flow on S^ = 1 implies that the a-limit of the orbit in 5i with 
Si 7 ^ —1 must be one of the transversally hyperbolic sources. 

This establishes Lemma 101 


Dynamics on KP 

The invariant subset KP is defined by setting z = 0 and S^ = 1 ; it can be represented by 
the Cartesian product of the circle (S^ = 1) in the Sj-space times the s^-triangle given by 
{0 < si, S 2 , S 3 < 1 , ~ 1 }' space possesses a simple structure: since 
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Figure 4: Schematic representation of the flow of the dynamical system on the invariant 
set KP = {z = Q)r\ = 1), which is the Cartesian product of the = l)-circle and the 
Sj-triangle. The depicted fixed points are the Kasner fixed points and the lines of fixed 
points TL^ (<-> Sj = 2) and QL^ (<-> Sj = —2). All orbits are heteroclinic, like the orbit 
through the arbitrary point P that connects KC 2 with KC 3 . 


S' = const for all orbits, the dynamical freedom resides in the Sj-spaces. 

A schematic depiction of the flow on KP is given in Figure 0 All fixed points are located 
on the boundaries of KP, i.e., on si = 0, S 2 = 0, or S 3 = 0. The vertices of the Sj-triangle 
are the Kasner circles KC°. If (Si,S 2 ,S 3 ) G (S^ = 1) is such that S^ = 2 (respectively 
Sfc = — 2 ) for some k, then the side Sfc = 0 of the triangle is a line of fixed points, TL^ 
(respectively QL®). Note that all fixed points are transversally hyperbolic on KP, and that 
they constitute the a- and cu-limit sets for all orbits on KP . The character of the fixed 
points, i.e., whether they are (transversal) attractors or repellors, depends on the sector 
of the circle (S^ = 1), see Figure 0] 

The results about the global dynamics on 5)^ and KP will turn out to be an integral part 
in the proof of Theorem uni which we will address next. First, we will prove the massless 
case of the theorem. 

Dynamics on 

Let 7 be an arbitrary orbit in the interior of Z ^, 7 7^ {F*^}. The function M( 2 ) is strictly 
monotonically decreasing on Z^ (except at F^, where it has a minimum), cf. (12011 ff.: hence 
we can use the monotonicity principle: the a-limit set 0(7) of 7 must be located on the 
boundaries of Z^^ i.e., on Sf or /C®. The first step in our analysis is to prove that the 
interior of the subsets 5? and JC^ cannot belong to 0(7), unless 7 is one of three special 
orbits. 
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Recall from our analysis of 5? that the fixed point G 5? is a hyperbolic sink on 5?. 
In the orthogonal direction, however, we obtain = 3. It follows that is a 

hyperbolic saddle in the state space and that there exists exactly one orbit 6 ^ that 
emanates from into the interior of Z^. iTheorem 14.21 implies that <5^ converges to the 
global sink as r —> oo.) 

Henceforth, let 7 be different from 5?. In order to show that a(j) does not contain any 
point of the interior of Sj^, we perform a proof by contradiction: assume that a(j} contains 
a point P of the interior of 5?; then the whole orbit through P and the tu-limit ti;(P) (as 
well as the a-limit) of that orbit must be contained in 0 ( 7 ). As already shown, is the 
global attractor on S^, hence Q!( 7 ) 3 a;(P) = D°. Since the saddle is in q;( 7 ), the unique 
orbit (5? emanating from it is contained in a(j} as well. Thus, ultimately, i.e., the 

point F®, must be contained in q;( 7 ); this is a contradiction, since F*^ is a sink. Therefore, 
7 cannot contain any a-limit point in the interior of 5?. We will now use similar reasoning 
repeatedly. 

Assume next that a(j) contains a point P in the interior of /C^, i.e., a point with = 1, 
0 < Si < 1 Vi. Suppose first that / ±2 for all k. Since P is an element of 0 ( 7 ), the 
whole orbit through P and the a-limit a(P) of that orbit must be contained in 0 ( 7 ). From 
the dynamics on /C^, cf. Figure it follows that a(P) is one of the Kasner fixed points 
on KC^, where i corresponds to the direction determined by Sj = max^Efc; we hence 
denote a(P) as Kp. Since 7 ^ ±2, it follows from the previous analysis that Kp is a 
transversally hyperbolic source on the subspace JC^] = 2 — ^^ WkZk > 0 yields 

that Kp is a transversally hyperbolic source on the whole space Z^. Since 0 ( 7 ) contains 
the transversally hyperbolic source Kp, that fixed point necessarily constitutes the entire 
a-limit set, i.e., 0 ( 7 ) = Kp. This is in contradiction to our assumption a( 7 ) 9 P. The 
omitted cases Sj = ±2 for some i will be dealt with next. 

Suppose that Sj = —2 for one index i. Assume that P lies in a( 7 ), therefore a(P) is 
contained in 0 ( 7 ) as well. The dynamics on JC^ implies that a(P) is a fixed point Qp 
on QL?, cf. Figure 01 This point is a transversally hyperbolic sonrce; = 1 in 

this case. By the same argument as above we obtain a contradiction to the assumption 
a( 7 ) 3 P. 

Finally suppose that Sj = 2 for one index i. When we assume that P is in 0 ( 7 ), then 
the cu-limit w(P) is contained in 0 ( 7 ). From Figure 0 we see that a;(P) is a fixed point 
Tp on TL?. The point Tp is a transversally hyperbolic saddle, since = 3, and 

there exists exactly one orbit that emanates from it, namely the orbit that connects Tp 
with Dj in 5?. Since Tp G 0 ( 7 ), that orbit must also be contained in a( 7 ). This is in 
contradiction to the previous result: a( 7 ) cannot contain interior points of 5^?. Hence our 
assumption 0 ( 7 ) 3 P was false: the a-limit of 7 cannot contain any interior point of JC^. 

Our analysis results in the following statement: There exist four special orbits, one trivial 
orbit corresponding to the fixed point F® and three orbits, the orbits (5®, that converge to 
the fixed points G 2?^. The a-limit set of every orbit 7 in Z^ different from F^ and 
(5? must be located on the boundaries of the spaces 5? and /C°, i.e., on the union of the 
boundaries of the cylinders Sf, which we denote by dS^ = dSf U dS^ U dS^. The set dS^ 
is depicted in Figured it comprises the lateral surfaces of the cylinders and the base/top 
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KC§ 



Figure 5: The heteroclinic network HP that exists on the set 55®. Its building blocks are 
the heteroclinic cycles W®, W®, W®. 


surfaces. 

All fixed points on 55® are transversally hyperbolic except for the points T®j: TL® consists 
of transversally hyperbolic saddles; in contrast, the fixed points on QL® are transversally 
hyperbolic sources; points on KC® with Sj > 1 , S, / 2 are sources while those with 
Sj < 1 are saddles. Combining the analysis of the preceding sections, see Figs. 01 and 
we obtain, more specifically: each point on QL® is a source for a one-parameter family of 
orbits that emanate into the interior of 2^®, and each point on KC® with Sj > 1 (Ej ^ 2 ) 
is the source for a two-parameter family. (The points with Sj = 1 on KC® are the two 
points Q®j G QL® and Q®^ G QL^- Each of these two points is a transversally hyperbolic 
source for a one-parameter family of orbits, however, those orbits are not interior orbits, 
but remain on the boundary of .) 

The non-transversally hyperbolic fixed points T®j are part of a special structure that is 
present on 55®: the set 55® exhibits a robust heteroclinic network ?f® (of depth 1 ), see, 
e.g., HH for a discussion of heteroclinic networks; the network ?f® is depicted in Figure 0 
a schematic depiction is given in Figure |H1 In particular we observe that the heteroclinic 
cycles ?f® of the spaces 5® appear as heteroclinic subcycles of the network. 

A straightforward analysis of the flow on 55® using the same type of reasoning as above 
leads to the result that there exist no other structures on 55® that could serve as a-limits 
for an interior 2^®-orbit 7 . We have thus proved the following statement: The a-limit 
of 7 is one of the transversally hyperbolic sources listed above, or it is the heteroclinic 
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network (or a heteroclinic subcycle thereof). This concludes the proof of the massless case 
of Theorem oi 

Global dynamics in the massive case 

Let 7 be an arbitrary orbit in the interior of the state space X. The function is 
strictly monotonically increasing on X (and on 1C), cf. ( 1181 ) ff.: moreover, vanishes for 
z ^ 1 and Sj ^ 0 (unless z —>■ 0 simultaneously). Hence, by applying the monotonicity 
principle we obtain that the a-limit set 0(7) of 7 must be located on including its 
boundaries. 

Consider the fixed point F® £ Z^. By Theorem 14.21 this fixed point is a global sink on 
Z^. In the orthogonal direction, however, we have z~^z'\-po = 2. It follows that is a 
hyperbolic saddle in the space X and that there exists exactly one orbit (j) that emanates 
from F*^ into the interior of X. ('Theorem 14.21 implies that 4> converges to FS^ as r ^ 00 ; 
thus, cj) represents the unique solution of the Einstein-Vlasov equations that isotropizes 
toward the past and the future.) 

Let 7 be different from </>. Assume that 0 ( 7 ) contains a point P of the interior of Z^; then 
the whole orbit through P and the w-limit ti;(P) must be contained in 0 ( 7 ). Theorem 14.21 
implies t(;(P) = F°, hence F° G «( 7 )- Since the saddle F*^ is in 0 ( 7 ), the unique orbit cj) 
emanating from it is contained in 0 ( 7 ) as well. Thus, ultimately, uj{4>), i.e., a point on 
FS^, must be contained in 0 ( 7 ); this is a contradiction, since FS^ consists of transversally 
hyperbolic sinks. We conclude that 7 cannot contain any a-limit point in the interior of 
Z^. 
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Since a( 7 ) must be located on the boundary on i.e., on Sf or the proof can be 
completed in close analogy to the proof in the massless case. We thus restrict ourselves here 
to giving some relations that establish that the sources on Z^ generalize to sources on X: 
on KC° we have z~^z'\y^qo = 2(1 + Sj), which is positive for all Sj > —1 and thus for Sj > 1 
in particular; for QL? we obtain z~^z'\qiQ = 4. We note further that z~^z'\y^o = 3; thus 
possesses a two-dimensional unstable manifold. (Orbits in that manifold converge to 
FS^.) Finally, note that along the heteroclinic cycle 7ii : T 22 —> T 32 —> T 33 —> T 23 —> T 22 , 
we obtain that z~^z' equals 6 ^ 2 , 2(1 -|- S 3 ), 2(1 -|- S 2 ), 6(1 — S 2 ), respectively; hence z~^z' 
is non-negative along the heteroclinic network Ti. 

This concludes the proof of Theorem 14.31 


5 Concluding remarks 

In this article we have analyzed the asymptotic behaviour of solutions of the Einstein- 
Vlasov equations with Bianchi type I symmetry. To that end we have reformulated the 
equations as a system of autonomous differential equations on a compact state space, 
which enabled us to employ powerful techniques from dynamical systems theory. 

Based on the global dynamical systems analysis we have identified all possible past at¬ 
tractors of orbits — both in the massless and massive case. We have found that an open 
set of solutions converges to the Kasner circle(s); in particular, for these solutions, the 
rescaled matter quantity satisfies 11 —> 0 toward the singularity, so that “matter does 
not matter.” However, we have seen that there exists an interesting structure that might 
complicate matters: there exists a heteroclinic network Ti.^ that might be part of the past 
attractor set. For solutions that converge to H has no limit toward the singularity, 
since H 7 ^ 0 along parts of the network i.e., matter does matter for such solutions. It 
is not clear at the present stage whether the set of solutions converging to is empty, 
or, if non-empty, of measure zero or not (the flow on the boundary subsets gives a hint 
that it might be a three-parameter set, i.e., a set of measure zero). In any case TL^ will be 
important for the intermediate dynamical behaviour of some models, and thus there are 
significant differences between Bianchi type I perfect fluid models and models with Vlasov 
matter. 

If a generic set of orbits converges to Ti.^, this will have considerable consequences. Bianchi 
type I perfect fluid models play a central role in understanding the singularity of more 
general spatially homogeneous models pQ, as well as general inhomogeneous models [HI. 
The importance of Bianchi type I perfect fluid models is due to Lie and source “con¬ 
tractions” in spatially homogeneous cosmology and the associated invariant subset and 
monotone function structure (which is quite similar to the hierarchical structure we have 
encountered in the present paper), and asymptotic silence in general inhomogeneous mod¬ 
els [ 2 ]. Similarly, we expect that Bianchi type I Einstein-Vlasov models hold an equally 
prominent place as regards singularities in more general — spatially homogeneous and in¬ 
homogeneous — Einstein-Vlasov models. Hence the resolution of the problem of whether 
the heteroclinic network attracts generic solutions determines if Einstein-Vlasov models 
are generically different from general relativistic perfect fluid models in the vicinity of a 
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generic singularity, or not. 

In this article we have not considered a cosmological constant, A. The effects of a positive 
cosmological constant can be outlined as follows: since p —> oo toward the singularity, it 
follows that A can be asymptotically neglected and hence that the singularity structure is 
qualitatively the same as for A = 0. However, toward the future A destabilizes FS^, which 
becomes a saddle, and instead solutions isotropize and asymptotically reach a de Sitter 
state. 

We conclude with some remarks on different formulations. We have seen that the variables 
we have used to reformulate the equations as a dynamical system yielded multiple repre¬ 
sentations of some structures, e.g., the Kasner circle. Replacing Si with Ei = /H, i.e., 

the Hubble-normalized spatial frame variables of PHH, and using y = m?H ^ instead of 
z, yields a single Kasner circle. The latter variables, however, are not bounded; indeed, 
they blow up toward the future in the present case. It is possible to replace the variables 
by bounded variables, however, variables of this type lead to differentiability difficulties 
toward the singularity. Issues like these made the variables we employed in this article 
more suitable for the kind of analysis we have performed. However, Fij-variables, or “Hj- 
based” variables would have been more suitable to relate the present results to a larger 
context; but it is not difficult to translate our results to the Fij-variables variables used 
in uma, where the relationship between the dynamics of inhomogeneous and spatially 
homogeneous models was investigated and exploited. 
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A Dynamical systems 

In this appendix we briefly recall some concepts from the theory of dynamical systems 
which we use in the article. 

Consider a dynamical system defined on an invariant set X C M™. The cu-limit set u){x) 
[a-limit set a(a:)] of a point x £ X is defined as the set of all accumulation points of the 
future [past] orbit of x. The simplest examples are fixed points and periodic orbits. 

The monotonicity principle ^ gives information about the global asymptotic behaviour 
of the dynamical system. If M : X ^ M is a function which is strictly decreasing along 
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orbits in X, then 


(28a) 

(28b) 


for all X £ X. 


ujix) c e x\x 

a{x) C E X\X 


lim M{C) / supM} 

X 

lim M(0 7 ^ inf M} 


Locally in the neighbourhood of a fixed point, the flow of the dynamical system is deter¬ 
mined by the stability features of the fixed point. If the fixed point is hyperbolic, i.e., if the 
linearization of the system at the fixed point is a matrix possessing eigenvalues with non¬ 
vanishing real parts, then the Hartman-Grobman theorem applies: in a neighbourhood 
of a hyperbolic fixed point the full nonlinear dynamical system and the linearized system 
are topologically equivalent. Non-hyperbolic fixed points are treated in centre manifold 
theory: the reduction theorem generalizes the Hartman-Grobman theorem; for further 
details see, e.g., m- If a fixed point is an element of a connected fixed point set (line, 
surface,...) and the number of eigenvalues with zero real parts is equal to the dimension 
of the fixed point set, then the fixed point is called transversally hyperbolic. Application 
of the centre manifold reduction theorem is particularly simple in this case. (The situation 
is analogous in the more general case when the fixed point is an element of an a priori 
known invariant set that coincides with the centre manifold of the fixed point.) 


B FRW and LRS^ symmetry 

In this section we discuss in detail the sets FRW and LRS,, connected with solutions 
exhibiting FRW or LRS geometry. 

To begin with, we prove that the fixed point F^ on z = 0 is well-defined. Since the 
defining equations for F° are wi = W 2 = W 3 = 1/3, we must show that these equations 
indeed possess a unique solution (si, S2, S3) for all distribution functions /q. Setting z = 0 
in m implies that equations wi = W 2 = VJ 3 = 1/3 are equivalent to the system 

« := y /o [si'i'i - S 2 V 2 ] Skvl^ ^ d^v = 0 (29) 

and u = 0, where v is defined by replacing [siuf — S 2 V 2 ] by [siv\ — S 3 U 3 ] in On the 

three boundaries of the space {(si, 52 , 53 ) | s* > = 1} the functions u and v are 

monotonic; their signs are given in Figured The derivative du/dsi is manifestly positive, 
du/ds 2 is negative, hence gradrt is linearly independent of the surface normal ( 1 ,1, 1 ), and 
it follows that u = const describes a curve for all const E M. The same argument applies 
to V, since dv/dsi > 0 and dv/dss < 0. Figure d reveals that n = 0 (u = 0) connects 
the upper (right) vertex of the ( 51 , 52 , S 3 )-space with the opposite side. Investigating 
(gradtt — A grad u) we find that the first component is manifestly positive when A < 2/3 
and negative when A > 3/2, the second component is negative when A < 3, and the third 
component is positive when A > 1/3, which implies that (gradu — A grad u) is linearly 
independent of the surface normal (1,1,1) for all A. It follows that all equipotential curves 
of the functions u and v intersect transversally; hence u = 0 and u = 0 possess a unique 
point of intersection, which proves the claim. 


21 


S3 



Figure 7: The functions u and v are monotonic along the boundaries of the space 
{(si, 52, S 3 ) I Si > 0, ■Sfc = !}• 


By establishing existence and uniqueness of the fixed point for all /o, we have shown 
that for all distribution functions there exists a unique FRW solution of the massless 
Einstein-Vlasov equations. 

The situation is different in the massive case. A FRW solution is characterized by the 
equations Sj = 0 Vi, wi = W2 = W3 = w, since this yields s, = const (and a rescaling of the 
spatial coordinates then results in gij oc 6 ij.) However, for a general distribution function 
/o, these equations are incompatible with the Einstein-Vlasov equations; in other words, 
the straight line Sj = 0 Vi, wi = W 2 = W 3 = w is not an orbit of the dynamical system. 
Hence, in the massive case, the Einstein-Vlasov equations do not admit a FRW solution 
for arbitrary /q; the distribution function /o is required to satisfy FRW compatibility 
conditions, see below, in order for a FRW solution to exist. 

Note, however, that for each /o, there exists exactly one orbit that originates from F*^ 
and ends on FS^, see Section^ i.e., there exists a unique solution of the Einstein-Vlasov 
equations that isotropizes toward the past and toward the future. This anisotropic solution 
can be regarded as a generalized FRW solution; if /o is compatible with the FRW geometry, 
then the generalized FRW solution reduces to an ordinary FRW solution. 

The treatment of the LRS case is analogous: the subset LRSi (and, analogously, LRS 2 , 3 ), 
defined through the equations T ,2 = W 2 = W 3 , describes solutions exhibiting LRS 
geometry. (For a solution on LRSi, Equation (HSEl) entails S 2 {t) oc S 3 (r); by rescaling 
the coordinates one can achieve §22 = 933 , he., a line element in an explicit LRS form.) 
However, for general /o, the set LRSi is not invariant under the flow of the dynamical 
system. Consequently, for general /o, the Einstein-Vlasov equations do not admit solutions 
with LRS geometry. 
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More specifically, consider 


(S 2 - S 3 )' = -Sn [i(l - a;)(S 2 - S 3 ) - {W 2 - tng)] . (30) 

Hence, (S 2 — S 3 )' vanishes when S 2 = S 3 and W 2 = W 3 . From (CSEl) and (I13cl) we obtain 

an equation for re', 

w[ = — 2 wi 

where we have defined 

(1 - z)^ f fo [z + il- z) Y.k Skvl] dvidV2dV3 

= - - — 

J fo[z + { 1 - z) Yjk ^kvf\ ' dvidV2dV3 

note that Wi = fdf'^. Equation m implies 

(^2 - u;3)' = - (Si - S 2 ) (/3™ - Pi°i) - ^ (Si + 1) (/3« - /3«) , (33) 

when S 2 = S 3 and W 2 = W 3 . We conclude that the set S 2 = S 3 , iC 2 = W 3 is an 

invariant set of the dynamical system, iff W 2 = W 3 implies = P 33 and P 2 '' — (^ 3 '’■ 

(In the massless case, only the first condition is required.) These conditions are violated 
for general distribution functions; for the condition to hold /q must be of a certain type 
that ensures compatibility with the LRS symmetry. This is the case, for instance, when 
there exist constants 02 > 0 , 03 > 0 , such that /o is invariant under the transformation 
V 2 ( 03 / 02 )^ 3 , V 3 {a 2 /a 3 )v 2 \ e.g., /o = h{vi,vlvl), or /o = /o(ui,a^i;^ + a^ul); in 
the latter case W 2 {t) = W 3 {t) implies a\s 2 {T) = a|s 3 (T). 

Note finally that a distribution function /o is compatible with a FRW geometry, if it is 
compatible with all LRS symmetries. This means, that for instance /o = /o(ai'i^i + 

a|i;|) is compatible with the FRW symmetry and thus admits a unique FRW solution of 
the Einstein-Vlasov equations. 


Si - ^ Sfc Qiufc + Vlfc ) + I {wi + /3^'>) 


(31) 


C Future asymptotics 

In this section we give the proof of Theorem 14.21 

Theorem HU The LO-limit of every orbit in the interior of the massive state spaee X 
[massless state space Z^] is one the fixed points ES^ [the fixed point F^[. 

Proof. Consider first the state space Z^ of massless particles and the associated sys¬ 
tem m- The function M( 2 ), cf. ( 1201 ) ff.. is well-defined and monotonically decreasing 
everywhere except for at the fixed point E*^, where it has a global minimum. On the 
boundaries 5? (given by s* = 0) and JC^ (S^ = 1) of the state space Z^, the function M( 2 ) 
is infinite. Therefore, application of the monotonicity principle yields that the w-limit of 
every orbit must be the fixed point F*^. 
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In the massive case consider (I13ap in the form 


S' = -3fi 


i(l - u))(l + Si) - ^(1 - 3w) - Wi 


(34) 


The r.h.s. is positive when Sj < —1 and z > 0 {w < 1/3). This implies that the hyperplanes 
Sj = — 1 constitute semipermeable membranes in the state space T, whereby the “triangle” 
(Si > —1) n (S 2 > —1) n (S 3 > —1) becomes a future invariant subset of the flow ([TTn) . 

The first part of the proof is to show that every orbit enters the triangle at some time Te 
(and consequently remains inside for all later times). 

Assume that there exists an orbit with Sj(T) < — 1 for all r (for some i). From (USEl) we 
infer that 

s' = —2sj [sj(Sj — Sj) + — Sfc)] > 0 (35) 

if Sj < —1, and that s' > 0 if Sj = —1; hence Sj(r) > Sj(To) = const > 0 for all r e [tq, 00 ). 
From m we obtain 


1 




= 1 - \wi{l + Sj) - \wj{l + T,j) - ^Wkil + Sfc) > 

—0 000 

> l — Wj — Wk = (1 — ^w) + Wi > const > 0 , 


(36) 


since Sj > const > 0. Consequently, n(r) > const > 0 for all r G [ro,oo). It follows 
from dni) that 

s' > const > 0 (37) 

for all r G [ro,oo) by the same argument. This is in contradiction to the assumption 
Sj < — 1 for all T. 

Thus, in the second part of the proof, we can consider an arbitrary orbit 7 and assume, 
without loss of generality, that 7 (r) lies in the S-triangle for all r G [re, 00 ). Equation (I13c|) 
leads to 

z' = 2z{l- z)Y,Sn{l + ^n)>0 (38) 

n 

for all T G [re,oo), hence z{t) > z(Te) > 0 for all r G [re,oo). 

We define the function N by 

iV = (I + Si)(l + S 2 )(l + S 3 ). (39) 


The derivative can be estimated by 


N' > 3ITV 



I - Zw 

1 + Sn 


(40) 


Since w{t) < const < 1/3 (because z{t) > const > 0), N' is positive when at least one 
of the Sj is sufficiently small, i.e., when N itself is small (a detailed analysis shows that 
N' > 3r2A'[—(3/2)(l —?/;) + V3{1 — 3t(;)iV“^/^]). We conclude that there exists a positive 
constant Nq such that N{t) > Nq for all r G [rg, oo). This in turn implies that there exists 
> 0 such that Sj(r) > —1-\-v for all i for all r G [re,oo), whereby z' > 2z{\ — z)v for 
all r G [re,oo). 
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It follows that the a;-limit of 7 must lie on z = 1, i.e., on . Taking into account the 
simple structure of the flow on Z^ ^ characterized by fl' = 3(1 — we conclude that the 
fixed points FS^ given by Si = S 2 = Sa = 0 are the only possible w-limits. ■ 

Remark. In order to demonstrate the versatility of the dynamical systems methods, we 
have chosen here to prove Theorem I4.2l bv using techniques that are slightly different from 
those employed in Section |1] (which exploit the monotonicity principle). However, it is 
straightforward (in fact, even simpler) to give a proof by making use of the hierarchy of 
monotone functions. Indeed, the function M(i) ensures that the w-limit of every orbit lies 
on Z^ or Si] modulo some subtleties, we can exclude that Si is attractive by using the 
monotone function and the local properties of the fixed points. 


D The spaces Sf — interpretation of solutions 

The flow on the boundary subsets 5? is of fundamental importance in the analysis of the 
global dynamics of the state space, see Section 14.21 Note that except for all attractors 
(D?, QL^, KC^, and the heteroclinic network) lie on Sf. For a depiction of the flow on 
5® see Figure |31 In the following we show that orbits on 5)^ represent solutions of the 
Einstein-Vlasov system that are associated with a special class of distribution functions. 
Furthermore, we investigate in detail solutions that converge to the subcycle of the 
heteroclinic network. 

Consider a distribution function /o of the form 


Mvi,V2 ,V3) = 5(ui)/o®'^(u2,U3) , (41) 

where fQ^^{v 2 ,V 3 ) is even in V 2 and ^ 3 . In the case of massless particles, m = 0 (and z = 0 
respectively), we obtain 

gJJ f yred ^2 ^^ 22.^2 ^33^2j “1/2 

w,=0, = ^- (j = 2,3), (42) 

/ dvidv2dv3 


where and can be replaced by S 2 and S 3 , if desired. In the unbounded variables g®* 
the equations read 


S; = -H[l + Si], = -2g^\l + Zi) (43a) 

S' =-H[l + S, - 3w;,] , {g^iy = -2g^\l + ^,) (j = 2,3), (43b) 

cf. the remark at the end of Section [2| In particular we note that the equation for g^^ 
decouples; hence the full dynamics is represented by a reduced system in the variables 
(Si, S 2 , S 3 ,( 7 ^^), which coincides with the system (IlHl) on the invariant subset g^^ = 0. 
In analogy to the definitions (t7^ we set 


Si = 0 , 


.,22 


S2 = 


g22 _|_ g33 


S3 = 


9 


33 


g22 g 


33 


(44) 
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so that S2 + S3 = 1 . This results in the dynamical system 

= —n[l + Si] , Si = 0 (45a) 

s' = - 0[1 + Sj - 3wj] , s' = -2sj[Sj - (S 2 S 2 + S 3 S 3 )] {j = 2,3) . (45b) 

This system (|i 5 |) coincides with the dynamical system m induced on 5 i (which is ob¬ 
tained by setting z = 0, thus w = 1/3, and si = 0 in (H^ i. 

Our considerations show that the flow on 5 i possesses a direct physical interpretation; 
orbits on 5 i represent solutions of the massless Einstein-Vlasov system of Bianchi type I 
with a “distributional” distribution function of the type dm. Note that the system dm 
on 5 i must be supplemented by the decoupled equations (lllbl) and {g^^Y = —25(^^(1-|-Si) 
in order to construct the actual solution from an orbit in 5 i. 

Two structures in 5 i are of particular interest: the fixed point D)* and the heteroclinic 
cycle Til, see Figure 01 The fixed point represents an LRS solution (associated with a 
distributional /o); it is straightforward to show that the metric is of the form 

gu = const , 522 oc , 533 oc , ( 46 ) 


and H = ( 4 / 9 )t 

The orbit T22 —> T32, which is part of corresponds to a solution 

511 = 5ii ) 522 = 522(3^^0*)^ ) 533 = 533 ; (47) 

here, H = ( 3 t)“^; Hq is a characteristic value of H. For the orbit T33 —> T23 the result is 
analogous with 522 and 533 interchanged. A more extensive computation shows that the 
orbit T32 —> T33 leads to 

5ii=5?i, 522 = 522 [log(l + 3i7oi)]^ , 533 = 533(1 + 37^0^)^ , ( 48 ) 

together with H = Hq {1 + 3HQt)~^{l + [log(l -|- 3HQt)]~^). (Note that 3Ht is always close 
to 1 and approaches 1 for t —> 0 and t ^ 00.) The result for the orbit T23 ^ T22 is 
analogous with 522 and 533 interchanged. 

Now consider an orbit converging to the heteroclinic cycle as r ^ —00, i.e., t \ 0 . Since 
the orbit alternates between episodes where it is close to one of the four heteroclinic orbits, 
we obtain a solution with alternating episodes of characteristic behaviour of the type dm 
and dm; transitions between the episodes correspond to the orbit being close to the fixed 
points. 

Let denote a monotone sequence of times such that the solution is in episode (n) 
at time (i.e., the orbit is close to one of the four heteroclinic orbits and far from 
the fixed points); \ 0 as n ^ 00. Since 3 Ht ^ 1 as t \ 0 , the sequence 
gives rise to a sequence defined by = 1 . During episode (n) the solution 

exhibits a characteristic behaviour of the type (14711 or (14811 with Hq = (and 522 = 522^ j 
533 ~ 533^)- ^ transition from one episode to another involves a matching of the constants. 

Example. Suppose that the orbit is close to the heteroclinic orbit T32 —> T33 in episode 
(n). We obtain a behaviour of the type (im with Hq = As gets small we see 
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that §22 ~ 933 ~ 533^ ■ The next (as t\ 0 ) episode corresponds to the orbit 

running close to T 22 —> T 32 ; the behaviour of the solution is (I47jl with and 

Hq = The transition between the episodes (n) and (n + 1) is thus straightforward; 

^' 22 '''^^-^^^"'''"^^)^ = and fl' 33 '''^^ = 533 ^ Matching episodes (n + 1) and (n + 2) 

is slightly more involved. The orbit is close to the heteroclinic orbit T 23 —> T 22 in episode 
(n + 2 ), where 


511 = 5ii , 


522=5S+'^(l + 3F(-+2)t)2, 


(n+2) 

933 = 933 


log{l + . (49) 


Close to transition time, when is large, we get §22 = t)"^ and 

933 = 533 ^^^(logThe transition between episode (n + 1) and (n + 2) thus 
involves that §33 begins to decay logarithmically from having been approximately constant. 
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